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We propose a new parameterization of the antineutrino flux from core coUapse supernovae, that 
aUows an interpretation of its astrophysical parameters within the Bethe and Wilson scenario for 
the explosion, and that leads to a reasonable (smooth) behavior of the average energy and of 
the luminosity curve. We apply it to analyze the events observed by Kamiokande-II, 1MB and 
Baksan detectors in correlation with SN1987A. For the first time, we consider in the same analysis 
all data characteristics: times, energies and angles of the observed events. We account for the 
presence of background and evaluate the impact of neutrino oscillations. The hypothesis that the 
initial luminous phase of emission (accretion) is absent can be rejected at the 2% significance level. 
Without the need to impose external priors in the likelihood analysis, the best-fit values of the 
astrophysical parameters are found to be in remarkable agreement with the expectations of the 
standard core-collapse scenario; in particular, the electron antineutrino-sphere radius is 16 km, the 
duration of the accretion phase is found to be 0.55 s, and the initial accreting mass is 0.22 Mq. 
Similarly the total energy emitted in neutrinos is 2.2 x 10^^ erg, again close to the expectations. 
The errors on the parameters are evaluated and found to be relatively large, consistently with the 
limited number of detected events; the two dimensional confidence regions, that demonstrate the 
main correlations between the parameters, are also given. 

PACS numbers: 97.60.Bw Supernovae; 26.30.Jk Weak interaction and neutrino induced processes; 95.55.Vj 
Neutrino detectors; 14.60.Pq Neutrino mass and mixing. 



I. INTRODUCTION 

We begin recalling the interest of core collapse super- 
novae, the status of their understanding, and the expec- 
tations for neutrino emission in the standard scenario. 
Next, we discuss in Sect. II Bl the motivations for an im- 
proved analysis of SN1987A observations in the context 
of the standard scenario. Finally, we offer an outline of 
the present investigation. 



A. Neutrino emission in core collapse supernovae 

Core collapse supernovae (SN) are astrophysical events 
in which all known forces interplay with each other in ex- 
treme physical conditions. An adequate modeling of the 
processes occurring during this event would be impor- 
tant to obtain information on the left-over compact star 
[H , on nucleosynthesis 0, H, Q , on the properties of the 
supernova remnant [5l| , and on the expected signals dur- 
ing the explosion; in particular, gravitational waves and 
neutrinos @, 0] ■ 

Because of the complexity of the problem, the model- 
ing of the physical processes is still in evolution, but it 
is generally accepted that the role of neutrinos is critical 
for the energy transport as first suggested in 0]. The 



collapse and the formation of a compact object, like a 
neutron star, have to pass through substantial neutrino 
emission, see, e.g., @. The details of how the explosion 
takes place and how the neutrinos are emitted are less 
clear and necessarily model dependent. In this work we 
focus on the only mechanism that has been studied in 
some detail: the neutrino-driven mechanism also known 
as Bethe and Wilson scenario [l0| or delayed scenario 
for the explosion. In the neutrino-driven mechanism, 
the explosion of the massive star receives crucial assis- 
tance from the energy deposition due to an initial, in- 
tense neutrino luminosity. Although the viability of this 
mechanism cannot be considered fully demonstrated at 
present^, recent theoretical results ll^ilSj encourage the 
opinion that the neutrino-driven mechanism works for 
certain core collapse SN. 

In the neutrino-driven mechanism, there are two main 
phases of neutrino emission: 

i) A thermal phase, called cooling, occurring when the 
proto-neutron star cools quietly. This phase involves 
most of the emitted neutrinos, 80-90% in energy. 

ii) A brief and very luminous neutrino emission, here 
termed accretion, that should involve a lower amount of 
neutrinos, 10-20% in energy. The accretion phase charac- 
terizes the neutrino-driven mechanism of the explosion, 
and it is expected to occur in the first stage of neutrino 
emission. In this phase, the matter is rapidly accreting 
over the proto-neutron star through the stalled super- 
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nova shock wave. The two most important processes of 
neutrino emission are 



P 



and 



P Ve 



(1) 



due to the abundant presence of nucleons and of quasi- 
thermal e^e~ plasma. These types of neutrinos {ve and 
v^) transfer to the star a small fraction of their energy, 
/ ~ 0.1, necessary to revive the stalled shock wave. 
See for a wide description of the phase of accretion, 
enriched by analytical arguments. 

The neutrinos from phases i) and it) can be observed in 
conventional supernova neutrino detectors [l5l | (namely, 
water Cherenkov and scintillators). In particular electron 
antineutrinos give signal mainly through inverse beta de- 
cay reaction on free protons: 



Ve P 



(2) 



Thus the existence of the accretion and cooling phases, 
generically expected in the neutrino-driven scenario, can 
be experimentally verified. 



B. What can we learn from SN1987A observations? 

SN1987A is the first and the only occasion at present 
to test the credibility of the various hypotheses on how a 
SN works. In fact, the events observed by 1MB [H, 
Kamiokande-II [13, [l^ and Baksan ^201 , represent a his- 
toric opportunity to investigate the physics of the col- 
lapse and of the explosion. A very extensive literature 
testifies the effort to extract information from these data 

43|, |44||. Usually, a specific characteristic of the SN1987A 




data is studied, most frequently the energy distribution. 
The energy and the time distributions are jointly con- 
sidered in a few analyses, but most often describing the 
neutrino emission with (overly) simple models, a proce- 
dure that is only partially justified by the limited amount 
of events collected. 

The next neutrino observation will be an extraordinary 
occasion to progress, but recall that SN are rare on hu- 
man time scale and it is not possible to reliably predict 
when the next one will happen. Thus, we should try our 
best using the only data that we have at our disposal, and 
in particular, we should attempt to address the question 
on whether there is a hint of accretion from SN1987A 
observations, as expected. In this respect, a point that 
deserves to be stressed is that all detectors observed a 
relatively large number of events in the first second of 
data taking, about 40 %: there are in fact 6 events in 
Kamiokande-II, 3 events in 1MB and 2 events in Baksan. 

A milestone for the point of discussion and more in 
general for SN1987A data analysis is the paper of Lamb 
and Loredo [s^ (LL in the following), where it is argued 
that the SN1987A observations can be used to claim for 
an evidence of the accretion phase. The LL paper, widely 



cited in theoretical and experimental reviews, is generally 
considered a useful application of refined statistical tech- 
niques and, in the present paper, we will provide the first 
independent verification of their results. However, we 
deem that, in view of the importance of their claim and in 
the light of various advances in neutrino physics (e.g., in 
oscillations, [13, H \M, 113 and [Hi [H, [SJ, [Hi, Hl ) , it is 
necessary to offer a critical discussion of the assumptions 
of the analysis by Lamb and Loredo. More specifically: 

1) the likelihood structure can be enhanced including a 
more accurate detection cross section [s^, the informa- 
tion on the directions of the events, a different treatment 
of the background [il ]: 

2) the theoretical model for neutrino emission, required 
for the data analysis, can be improved and approached 
as much as possible to the real signal expected from the 
numerical simulations of neutrino emission. In this re- 
spect ja specific criticism was raised by Raffelt and Mi- 
rizzi |42| . who emphasized that the two-phase parame- 
terization used in LL analysis leads to a sudden jump 
of the average neutrino energy while passing from accre- 
tion to cooling. This behavior is rather different from the 
expectations of the numerical simulations. 

In short, our main tasks are to present a somewhat dif- 
ferent likelihood, to propose an improved parameteriza- 
tion for neutrino flux, to include oscillations, and finally 
to evaluate the impact of the various new points for the 
analysis of SN1987A observations. 



C. Layout of the work 

The structure of this paper is the following: in Sect. 
im we discuss the construction of the likelihood function 
involved in the analysis of the data set, underlining the 
improvements carried in the description of detection rate; 
in Sect. IIIII we propose a new parameterization for neu- 
trino emission, building it step by step and discussing 
its features, leaving the description of certain technical 
details to the appendix; finally in Sect. IIVI we draw the 
results of our analysis. 



II. LIKELIHOOD CONSTRUCTION 

In this section we describe the likelihood that we 
adopted to compare the observed events and the assumed 
flux for the neutrino emission, stressing the novelties and 
the technical improvements. 



A. Signal rate 

The signal in each detector, triply differential in time, 
positron energy Ef, and cosine of the angle 9 between the 
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antineutrino and the positron is: 
' d cos 9 



R{t, Ee, cos 9) ^ Np-^{E^, cos 6) $p^(t,^,)x 

dE, (3) 

X^d (cos 61) f]d{Ee) -j^, 



where Np is the number of targets (=free protons) in 
the detectors, crp^p is the inverse beta decay cross section 
(Eq. [2) , rid the-detector dependent-average detection ef- 
ficiency, is the angular bias —1 for Kamiokande-II and 
Baksan whereas for 1MB ^d{cosd) = 1 + 0.1 cos^ [13, 
and, finally, is the electron antineutrino flux, differ- 
ential in the antineutrino energy E^, and discussed later 
in this work. 

The expected number of signal events fis is the cru- 
cial ingredient, along with the expected number of back- 
ground events /Xh, to construct the Poisson likelihood 
/i" exp(— /i)/n!, where /i = iJs + fJ-b and where n is 
the number of observed events. To evaluate fis we use 
Eq. [31 the number of expected signals in a bin is just 
R{t, Ee, cos 0)dtdEed cos 9; the total number of the events 
is the integral of R{t, E^, cos 9) over its three variables. 

As mentioned in the introduction, our analysis is sim- 
ilar to the one of Lamb and Loredo .'!(> , with whom we 
agree within errors when we strictly stick to their pro- 
cedure. The signal rate that we adopted in this paper 
departs from their one in the following points: 



1. Cross section and event direction 

We adopt the inverse beta decay cross section calcu- 
lated in [56| and in particular we use the differential 
expression day^p/ dcos9 given in Eq. (20) of that pa- 
per. The energy of antineutrino is given in terms of the 
positron energy Ef, and of the angle 9 between the an- 
tineutrino and the positron directions: 



Ey = 



Ee+S. 



1 - {Ee -Pe COS 9) I rup' 



(4) 



where 5_ = {rn^^ — rrip — ml)/ {2mp) = 1.294 MeV andpe 
is the positron momentum. The new total cross section 
agrees at 10 MeV with the one used by Lamb and Loredo, 
whereas at 20 MeV (30 MeV) it is 6% (12%) smaller. 



2. Ejficiency 

Following the traditional approach and differently 
from [36] we include in Eq. [3]the detection efficiency as a 
function of the true energy of the event. A formal justifi- 
cation of our procedure is given in [4^ , that is in contrast 
with the formal justification in Appendix A of [36l |. 

Our procedure simply accounts for the evident fact 
that the expected number of signal events /Zc, should in- 
clude all relevant detector dependent features: loss of 



events due to light attenuation, fiuctuations of the num- 
ber of photoelectrons, detector geometry, etc.. These fea- 
tures produce an imperfect {rjd < 100%) detection effi- 
ciency, that means that only a fraction of the produced 
positrons is actually detected. We have in mind an av- 
erage efficiency evaluated by a MC procedure, namely 
1) simulating several events with true energy Ee but lo- 
cated in the various positions and emitted in all possible 
directions, then 2) counting the fraction of times that an 
event is recorded, finally 3) deducing also the 'smearing' 
(=average error as a function of Eg). 

It is possible to argue in favor of our procedure by 
considering the following situation. Imagine two detec- 
tors that differ by the detector efficiency: rj = 100% in 
the first one, rj = 10% in the second, and with a rate of 
positron production equal to the background rate. Sup- 
pose that each of these detectors observed one event. Ac- 
cording to the procedure of [36*1 (used, e.g., to obtain their 
Tab. VI) the probability that the observed event is due 
to a signal is 50% in both detectors, that we find para- 
doxical: the better the detector, the higher should be the 
chances of a signal. Instead, adopting our procedure, the 
probability that the event is due to a signal is 50% in the 
first detector and 9% in the second one, which we find 
more plausible. 

We recall (in agreement with [S^H^I) that for an even 
more refined analysis of the data, one should not use the 
average detection efficiency, but should rather evaluate 
the specific detection efficiency and background rate for 
any individual event. In our understanding, a correction 
on individual basis of this type was performed only to 
assess the errors on the energies of the events, see [191] . 



B. The assumed likelihood 

We estimate the theoretical parameters by the x^: 

X^^~2 J2 log('C<i), (5) 

d—k,i,b 

where Cd is the likelihood of any detector {k,i,b are 
shorthands for Kamiokande-II, 1MB, Baksan). We use 
Poisson statistics; dropping constant (irrelevant) factors, 
the 'unbinned' likelihood of each of the 3 detectors is: 



X [f + jR{t,, Ee,cos9,)C^{Ee)dEe 



(6) 



We denote by R{t) the integral of R{t, Ee, cos 9) over the 
variables Ee and cos6'. In 1MB, the live-time fraction is 
fd = 0.9055 and the dead-time is Td — 0.035 s, whereas 
for the other detectors fd = ^ and Td ~ 0. Each detec- 
tor saw Nd events; their time, energy and cosine with 
supernova direction are called ti, Ei and Ci (i — l...Nd). 

Ci is a Gaussian distribution that includes the esti- 
mated values of the energy Ei and the error of the en- 
ergy 5Ei for each individual event, accounting for the 
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detector-dependent effects on energy measurements. The 
inclusion of the error on the measurement of cos 9 does 
not change significantly the likelihood, so we simply set 
COB 9 = Ci for each event. ^ Finally, we do not include an 
error on the event times 5ti, since the relative time of 
each event is precisely measured. 

The absolute times have not been measured precisely 
enough (except for 1MB). However, the procedure of 
analysis that we adopt requires only the relative times 
between the events: the experimental input is 5ti = 
f^^ — . The times ti are defined as follows: 



(7) 



where > is the offset (or delay) time between the 
first neutrino that reached the Earth (that, by definition, 
occurred at i = 0) and the first event that was detected. 
We introduce one parameter t°" for each detector, and 
fit their values from the data. The integral over the time 
in the first exponential factor of Eq. [5] is performed from 
the moment when the first neutrino reaches the Earth till 
the end of data taking, t = 30 s; the condition that all 
the data are included imposes mild restrictions, such as 
^Kii < 6 s, that do not have a relevant role in the analysis. 



1. Background 

The probability that an event is due to background is 
denoted by Bi in Eq. [6l It is calculated as Bi = B{Ei): 
Bi is the measured background rate for the given en- 
ergy [Hz/MeV]. The background distribution differential 
in time, energy and cosine is B{Ee)/2; the factor 1/2 
describes a uniform cosine distribution. This definition 
is different from the one of LL, Bi = jB{Ee}Ci{Ee)dEe, 
that has been argued to be inaccurate in j4J|. The val- 
ues of Bi that we use for Kamiokande-II are given in 
Appendix A of . The events of Kamiokande below 7 
MeV have a higher background rate than found by LL, 
those above 9 MeV a lower background rate, while the 
other ones stay almost unchanged. The changes for Bak- 
san are instead negligible. It is fair to assume in good 
approximation that, in the time window of interest, 1MB 
observations are safe against background contamination. 



are needed to quantify the duration of the emission pro- 
cess, the intensity of the emission and the average energy 
of the antineutrinos. The adopted time distributions are 
constructed to enforce the continuity of the instantaneous 
luminosity and of the average energy as found in the nu- 
merical simulations. Moreover, we tried to maintain the 
parameterization as simple as possible. 



A. Parameterized antineutrino fluxes 

1. Cooling phase 

In the last phase of the SN collapse the nascent proto- 
neutron star evolves in a hot neutron star (with radius 
Rns) and this process is characterized by a neutrinos and 
antineutrinos fiux of all species. This is the cooling phase; 
we use the suffix c in the corresponding symbols. 

A rather conventional parameterization of the electron 
antineutrino fiux, differential in the energy is: 



1 



[At:RI g^SEu.T.it))] (8) 



where the Fermi-Dirac spectrum of the antineutrinos is 

El 



g,AE.,T,{t)) 



(9) 



l + cxp[S,/Tc(t)] 
The time scale of the process is included in the function 



r,(i) = rcexpH/(4T,)]. 



(10) 



Eq. [5] describes an isotropic emission of antineutrinos 
from a distance -D(= 50 kpc in the case of SN1987A). 
[We use the symbol $° rather than to emphasize that 
flavor oscillations have not been included yet]. 

The astrophysical free parameters are Rc, T^^ and 
Tc namely: the radius of the emitting region (neutrino 
sphere), the initial temperature, and the time constant 
of the process. We recall which are the generic expecta- 
tions: Rc ~ Rns = 10 - 20 km, Tc = 3 - 6 MeV, and 
Tc =few-many seconds. Rather than using these a priori^ 
we will deduce the value of these parameters by fitting 
the SN1987A data, and later, we will compare the results 
with the expectations. 



III. ANTINEUTRINO FLUX DESCRIPTION 

In this section we describe the parameterization of the 
neutrino flux that we adopt. We describe the signal in- 
troducing three 'microscopic' (i.e., physically meaning- 
ful) parameters for each phase, that, roughly speaking. 



^ For the first 12 Kamiokande-II and for the 8 1MB events the value 
of cos S = Ci is measured; we set instead Ci = for the 5 events 
of Baksan and the last 4 events (out of 16) of Kamiokande-II. 



2. Accretion phase 

After the bounce, the simulations indicate that the 
shock wave, propagating into the outer core of the star, 
looses energy and eventually gets stalled. It forms an ac- 
creting shock that encloses a region of dissociated mat- 
ter and hot e^e~ plasma, where the weak reactions of 
Eq. [T] give rise to intense v^. and Ve luminosities. This 
emission lasts a fraction of a second. In the Appendix 
we describe in more details the conceptual scheme for 
emission: a neutron target exposed to a flux of thermal 
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positrons. The neutrons are treated as a transparent tar- 
get, for only a small fraction of antineutrinos is expected 
to couple with the star. This is the accretion phase; we 
use the suffix a in the corresponding symbols. 
The parameterized D,,, flux is 



Sttc 



{E,),Ta{t)) 



(11) 



where Nn(t) is the number of target neutrons assumed 
to be at rest and the thermal flux of positrons: 



ge+iE,+ ,Ta{t)) 



E^ 



1 + exp [E,^/Tait)] 



(12) 



is calculated at an average positron energy, namely 
Ee+{Eu) = \^Ej^f^ - In the energy range of inter- 
est, Eu = 5 — 40 MeV, a simple numerical approximation 
of the cross section for positron interactions is 



4.8 X IQ-^^El 
l + E^/{2m MeV) 



(13) 



The derivation of Eqs. llll and ll3l is given in the appendix. 

The average energy of the antineutrinos is roughly 
given by 5 Ta and the spectrum is slightly non-thermal,'^ 
mostly due to the presence of the cross section fie+n- For 
example, when Ta = 1.5, 2.5, 3.5 MeV, we get {E^) /Ta = 
5.5, 5.2, 5.0 resp ectively and 5 E^/Ta = 0.39, 0.41, 0.41 
where 5Ei, = ^ (El) — (E^)^. Another manifestation 
that the distribution is non-thermal is the scaling of the 
luminosity with the temperature, roughly as T^-different 
from the thermal scaling T^. In order to give some feel- 
ing of the antineutrino emission, if we suppose that at 
t = we have Ta = 2.5 MeV and Ma = 0.15 Mq (see 
just below for a precise definition of these quantities) we 
get a luminosity of 1.1 x 10^'^ erg/s; the same luminos- 
ity and average energy would be given by a black body 
distribution with Tc = 4.1 MeV and i?c = 82 km. 

There are two time dependent quantities in Eg. 1111 the 
number of neutrons Nn(t) and the positron temperature 
Ta{t). It is straightforward to introduce a temperature 
that interpolates from an initial value to a final value: 



Ta{t)=T, + iTf-T,) 



with 



T^Ta 

Tf = 0.6 Tc 



(14) 



where Ta denotes the positron temperature at the be- 
ginning of accretion (to be contrasted with T^, the an- 
tineutrino temperature at the beginning of the cooling 
phase). With this parametrization, the positron temper- 
ature reaches 0.6 at t — r„ , that is what is needed 



^ The deviation from the thermal distribution can be described 
by a 'pinching factor'— an effective chemical potential introduced 
to distort the Fermi-Dirac thermal spectrum-in the range 4-5, 
that decreases when Ta increases; e.g., for &Ei,/Ta = 0.41 the 
pinching factor is 4.2. 



match the average antineutrino energies, namely, to en- 
sure a continuous behavior of the average antineutrino 
energy (in particular at the end of the accretion phase 
and at the beginning of the cooling phase). The power 
m = 1 — 2 mimics the behavior found in numerical sim- 
ulations; we adopt m = 2 as a default value. 

Now we discuss the time evolution of the number of 
neutrons exposed to positrons N„{t), proportional to the 
luminosity in accretion. Our goal would be a luminosity 
that, at least for t ^ 0, decreases as 1/(1-1-^/0.5 s). This 
behavior is suggested by the numerical simulations and 
is advocated by LL [36[. However, when we allow the 
temperature to vary, we vary also the luminosity, that 
scales as NnT'^. Thus, we need to include an explicit 
factor {Ta{t) /Taf' ■ We arrive at our prescription for the 
number of neutrons exposed to thermal positrons: 



Nn{t) = ^ X Ma X 



jk{t) 



1 -l-i/0.5 s' 



(15) 



the fraction of neutrons being set to Yn = 0.6. Ma is the 
initial accreting mass exposed to the positrons thermal 
flux. The time-dependent factor 



jfc(t) =exp[-(t/T,)'=], 



(16) 



is included to terminate the accretion phase at i ~ r^. 
LL use fc = 10, which however leads to a very sharp drop 
of the luminosity at i ~ Tq. In our calculations, we will 
set instead fc = 2, a choice that offers the advantage of 
leading to a smooth (reasonable, continuous) luminosity 
curve, closer to the type of curves found in numerical 
simulations. [We will show in the next section luminosity 
curves with fc = 10 and with fc — 2.] 

Ultimately, the accretion phase involves 3 free param- 
eters: Ma, Ta and Tq, the same number of parameters 
of the cooling phase. Finally, we give rather generic ex- 
pectations on values of these parameters: Ma is certainly 
lower than the whole outer core mass (about 0.6 Mq); Ta 
is expected to sit in the few MeV range; and, finally, the 
accretion should last a fraction of a second, Tq ^ 0.5 s 
being a typical number. 



B. Temporal shift 

In the model advocated by Lamb and Loredo and 
adopted for the analysis of SN1987A data, the accretion 
and the cooling phases are contemporaneous. This has 
the consequence that for times t < Ta, the antineutrino 
distribution is not a thermal spectrum, but a composition 
(=sum) of two thermal spectra, whose average energies 
differ by more than a factor of 2 in the best-fit point. 
At low energies, the spectrum is dominated by the an- 
tineutrinos from accretion; at the highest energies, by 
the antineutrinos from cooling. The possibility to have 
a composite spectrum implies, among the other things, 
that it is easy (trivial) to reconcile the low energy events 
observed by Kamiokande-II and the high energy events 



6 



observed by 1MB in the first second: they simply be- 
long to different and simultaneous phases of emission. 
We will verify this statement later by a straightforward 
calculation-see Tab. HT] We note in passing that such a 
composite spectrum is used in the work of Lorcdo and 
Lamb [36] but this characteristic feature is neither com- 
mented or discussed there. The compositeness of the LL 
spectrum (i.e., the non-thermal tail) can be better per- 
ceived plotting it in logarithmic scale and/or by consid- 
ering the time-integrated spectrum, see Fig. 2 of [i^ . 

However, at the best of our knowledge there is no evi- 
dence from numerical simulations of a composite behav- 
ior of the spectrum; the instantaneous De spectrum is 
found to be quasi-thermal at any time and typically this 
property is shared also by the time-integrated spectrum; 
see, e.g., the discussions in [13, IH,!!^]. Deviations from 
a thermal distribution are observed, especially during ac- 
cretion [H, [iO, llll ; they can be effectively described by 
a 'pinching' parameter of the order of a few, that means 
that the high-energy tail of the spectrum is depleted-not 
enhanced as for the composite model advocated by LL. 

In short, we believe that, in absence of an explicit indi- 
cation from numerical simulations, the model that should 
be adopted in data analyses (and/or the null- hypothesis 
that should be tested) is the simplest one compatible with 
the numerical simulations, namely: a v,, spectrum quasi- 
thermal at any time-rather than, e.g., 'composite' [s^ . 
bimodal ^43| or exponentially decreasing in the energy 
[4^ . For this reason, we parameterize the antineutrino 
flux as follows: 



C. Neutrino oscillations 

Here we discuss the effects of neutrino oscillations on 
the observed v,. fluxes. The survival probability P for 
Ve emitted by supernova is dictated by two different in- 
teractions with the star medium. The first is the usual 
matter effect [13, [6l| of charged current interactions be- 
tween and the electrons of the matter. The second 
is the effect oi v — v interactions [H, [H, [13, Issl [6^ . 
that is known to be important in specific cases |67l . |68| 
and whose behavior has been quantified in certain ap- 
proximations [H, [zO, [zE Izll- In order to describe os- 
cillations we have to distinguish the two arrangements of 
the neutrino mass spectrum that are compatible with the 
present knowledge of neutrino properties (see, e.g., [TSj): 
For normal mass hierarchy the survival probability and 
the observed fiux are: 



= P 



(18) 



where we recall that is the fiux in absence of oscilla- 
tions. We have assumed that = $? and each term is 
the sum of the cooling fiux and accretion fiux. The v — v 
interaction is most relevant for inverted mass hierar- 
chy [ei, H [lE III, [Zl, [li- As an example the electron 
antineutrino survival probability given in 72] is: 



P= 1/2^(1 -P/) + U-3^/ 
= P <^l^ 



(19) 



$P,(t) = $a(i) + (1 - Jkit)) X $,(t - Ta). 



(17) 



where $a (resp., <I>c) is given in Eq. [TT] (resp., Eq. [8|) 
in the case when oscillations are absent. Recalling that 
jk{t) appears explicitly in Eq. [151 i-C-, in the accretion 
flux, it is clear that the effect of this function is simply to 
interpolate between the two phases of neutrino emission; 
in other words, the cooling phase begins around i = Tq. 

We note here a limit in which the likelihood becomes 
unphysical. Consider the case when only the first event 
detected by Kamiokande-II falls in the accretion phase. 
If the accretion were to last a very short amount of time, 
Ta ^ St2, and if t°^^ <C 6t2 — 0.107 s, the antineutrino 
fiux at the time of the first event t = 5ti — could be- 
come very large even if the number of expected number 
of events remains small. In this limit, the interaction 
rate in Eq. [3] becomes large, too. Thus, due to the fac- 
tor R{ti, Ee,cos6i) in Eq. [51 the likelihood can be made 
arbitrarily large, because the exponential factor in the 
same equation (that depends on the expected number 
of events) will not change much. The way to avoid this 
pitfall is to require a lower limit on Ta in the numerical 
calculations. A limit that is adequate for the analysis of 
data from SN1987A is > 0.3 s. We will demonstrate 
explicitly the presence of this unphysical behavior of the 
likelihood function later, when discussing its maxima. 



We adopt the usual decomposition of the mixing ele- 
ments in terms of the mixing angles: Ues = sin^ia and 
Uei = cos 6*12 cos 6*13 with 6*12 = 35° ± 4° and 6*13 < 10° 
at 99 % C.L. For the measured solar oscillation param- 
eters, the Earth matter effect is known to be small (see, 
c-g-, [ill)- We include it in the analysis anyway evaluat- 
ing the survival probabilities with the PREM model [t^] • 
In the case of normal mass hierarchy the probability 
P ^ 0.7 is reliably predicted and rather precisely known. 
Instead, for inverted mass hierarchy, P depends strongly 
on the unknown mixing angle ^13. In fact, the fiip prob- 
ability Pf (that quantifies the loss of adiabaticity at the 
'resonance' related to the atmospheric Am^) is: 



P/(S„0i3) = exp 



3.5 X 10- 



20 MeV^/^ 



(20) 

where the numerical value corresponds to the supernova 
profile Ne ^ 1/r^ given in [6l|. The predictions for the 
De fiux in the detector could be tested with a relatively 
large amount of data: a galactic supernova could turn 
out to be useful to discriminate the right mass hierarchy. 

We assume that an equal amount of energy goes in each 
species (equipartition hypothesis) during cooling. Fur- 
thermore, we suppose that the temperature of D^ and _Dj 
is in a fixed ratio with the D^ temperature. Following [78| 
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we assume: 

T{vr)lT{ve) = T{v^)lT{v,) = 1.2, (21) 

We tested that a value in the 1.0-1.5 or a deviation of 
the amount of energy stored in non-electronic neutrino 
species by a factor of 2 does not affect crucially the fitted 
antineutrino flux. In the accretion phase, we will suppose 
that only v,, and Ve are emitted in equal amount, whereas 
*^a(^a:) = ^he fit provides a reasonably fair description 
of the antineutrino flux anyway, but the estimation of 
the amount of energy emitted during accretion should be 
considered as a lower bound. 



IV. RESULTS 

In the previous two sections we constructed, step by 
step, a likelihood function that represents the probability 
function of the overall data set. This probability varies in 
the parameters space and depends on the model for the 
antineutrino emission. As anticipated, the main goal of 
this paper is to verify if the best-fit values of the parame- 
ters, that maximize the likelihood function, are physically 
acceptable. 

We remark that LL adopted a Bayesian analysis while 
we use a frequentist approach, which occasionally leads 
to some difference in the confidence levels, but not in the 
best fit points. More specifically to calculate the error 
on one parameter we use as a rule the profile likelihood, 
namely the likelihood evaluated fixing the parameter of 
interest and maximizing the other (nuisance) parameters. 
A similar procedure allows us to calculate the confidence 
regions. 

The structure of this section is the following: in 
Sect. IIV Al we describe the simpler (one component) 
model, in Sect. lTVB] we describe the more complete (two 
components) model. This second part illustrates the im- 
pact of each improvement in the description of the flux; 
the final result is discussed in detail in Sect. IIV B 51 



A. One component model 

Here we list the results of the likelihood maximiza- 
tion procedure when we include only the cooling phase. 
This model has 6 parameters: the temperature Tc, 
the duration Tc, the neutrinosphere radius Rc and the 
three detector offset times t°^^. It is identical to the "Ex- 
ponential cooling model" reported in [36l|, also termed 
"minimum" or "standard" model in [27[. The compari- 
son of the best-fit values of LL (i?c = 40 km, Tc = 3.81 
MeV, Tc = 4.37 s and t"^' = s (Hi) with our results 
{Rc = 44 km, Tc = 3.68 MeV, Tc = 4.43 s and f'' = s) 
is satisfactory, the agreement being at the level of 10%. 

We proceed and quantify the impact of the improve- 
ments in the likelihood described in Sect. [Ill 
The inclusion of the detection efficiency ri{Ec+ ) gives the 




2 2.5 3 3.5 4 4.5 5 5.5 6 

T,(MeV) 

FIG. 1: Two dimensional confidence regions for the cooling 
parameters Rc and Tc. The gray (red) contours are the expo- 
nential cooling result (68% and 90% C.L.) following Loredo 
and Lamb analysis. The dark contours are our result for one 
component cooling model. The best fit points are also dis- 
played. 

single most important effect. In fact we obtain: 

Rc = 30 km, Tc = 4.21 MeV, = 3.88 s. (22) 

Using the new cross section ap^p{E^^) for the inverse beta 
decay we get the new best-fit point values: 

Rc = 26 km, Tc = 4.58 MeV, Tc = 3.72 s. (23) 

Our assumption on the background Bi has a little effect 
on these values that become Rc = 26 km, Tc = 4.59 MeV 
and Tc — 3.81 s. 

Similarly the inclusion of the event directions, cos^^ 
(Eq. (20) in Isil), that produces: 

Rc = 26 km, Tc = 4.47 MeV, = 3.88 s. (24) 

The values of t°" parameters are always zero in these 
models, and it is easy to convince oneself that this is due 
to the fact that, in this one component model, the signal 
is forced to decrease with the time. 

The larger change concerns the radius Rc that dimin- 
ishes by 35% in comparison to the value given by Lamb 
and Loredo, approaching the expected neutron star ra- 
dius, but still twice larger than a typical value. In Fig.[T] 
we show the contour plots in the Rc — Tc plane, the 
two parameters that show the largest correlation among 
them. In gray we draw the 68% and 90% contours level 
that we obtain adopting the same likelihood function of 
LL, in black the contours level when we construct the 
likelihood function following Sec. |TT1 i.e., including all 
the structural improvements discussed above. 
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With the best-fit points it is possible to estimate the 
total energy emitted by neutrinos in this model. Hypoth- 
esizing that in the cooling phase all types of neutrinos are 
emitted, each one carrying away an equal amount of en- 
ergy, the total energy is: 

^.3,39xl0-f (25, 

This value should be comparable with the gravitational 
binding energy of the new born neutron star. Eh = (1—5)- 
Iq53 ^j^g result for exponential cooling model is 
El, = 5.02-10^'^ erg, namely a binding energy at the upper 
limit of this range, whereas our results E^ = 3.55-10^^ erg 
is included in the expected range. 



Sect. 


Rc 


Tc 




Ma 


Ta 


Ta 


t°"' 


Ax' 




[km] 


[MeV] 


[s] 


[Me] 


[MeV] 


[s] 


[s] 




El 


12 


5.46 


4.25 


5.59 


1.52 


0.72 


0. 


14.7 




14 


4.99 


4.76 


0.82 


1.75 


0.67 


0. 


11.2 




















IlIIBl 


14 


4.88 


4.72 


0.14 


2.37 


0.58 


0. 


9.8 


[mcl 


16 


4.62 


4.65 


0.22 


2.35 


0.55 


0. 


9.8 



TABLE I: The best-fit values of the astrophysical parameters 
for two components model neutrino emission. Each line of 
this table is an incremental step toward the final improved pa- 
rameterization. The last column shows the difference between 
the x' of our one-component (cooling) model and the x^ of 
each two-component model. 



B. Two components model 

Now we include both emission phases. Following the 
order of Sect. Ill, we improve, step by step, the emission 
model of the accretion phase. To describe accretion we 
add 3 new physical parameters, namely: the positrons 
temperature T^, the duration of the accretion phase 
and the initial accreting mass Ma- So the likelihood is a 
function of 9 parameters. The best fit results are given 
in table U and will be commented in details below. 

The last column of Tab. [J shows the values Ax^ = 
Xc ~ Xmi where is calculated with the best one- 
component model of the previous section, i.e., the last 
case of previous subsection (eq.[21]), and is calculated 
with the model described in the section indicated in each 
line of Tab. HI The larger the difference, the larger the 
evidence for accretion; a quantitative evaluation of the 
evidence, taking into account the increased number of 
parameters, is provided later. 

1. Effect of the new likelihood-improvements of Sect. [771 

The first line of Tab. U shows the best-fit results ob- 
tained using the likelihood function constructed in our 
Sect, [n] and one of the emission model used by LL, 
namely the model called "Exponential cooling and trun- 
cated accretion" (later called ECTA) [36j]. We see that 
the best-fit value for the initial accreting mass Ma is very 
large and hardly acceptable on physical basis: this pa- 
rameter is restricted by Ma < O.GMq for the reasons 
mentioned in Sect. IIII A 21 This result, however, is in 
agreement with what found by Lamb and Loredo, who 
fixed Ma = 0.5 Mq in all models of their table V, namely 
they inserted a "prior" in their analysis. This assumption 
reduces to 8 the number of free parameters and produces 
these best-fit values: 

i?c = 12 km, Tc = 5.40 MeV, Tc = 4.40 s, , . 
Ma = 0.5 Mq, Ta = 2.02 MeV, Ta = 0.70 s 

with Ax^ = 13.4. We call this best-fit point LL* to 
distinguish it by the true global maximum of likelihood 



function, shown in Tab. [H In passing, we note that our 
modifications of the likelihood do not really change the 
conclusions of Lamb and Loredo 179] . 



2. New spectrum and Ta{t) -improvements of Sect. \IIIA\ 

The second line of table H] shows that there is less 
need of the a priori on Ma when we exploit the correct 
parametrization of the accretion flux (accounting for the 
right kinematics of e"*"n process, using the new cross sec- 
tion and allowing the positron temperature to increase 
with the time). In fact the best fit value for the ini- 
tial accreting mass decreases to Mq = 0.82 Mq in this 
model, a value that is a bit larger than the reference 
outer core mass, 0.6 Mq, but now closer to the expected 
range. The accretion phase is characterized by a low 
mean energy and by a duration shorter than a second, 
as expected. The total energy carried in each phase is: 
Ec = 2.0 • 10^3 erg and Ea = 5.7 • lO'^^ erg, and the total 
binding energy is E^ = Ea + E^ — 2.5 • lO^'^ erg. It is 
important to note that, in this model and as in the LL 
model (see Sect. IIIIB|) the two emission phases are con- 
temporaneous. This implies that the emission spectrum 
is 'composite', a features that makes it easier to account 
for the difference between the average energies of 1MB 
and of Kamiokande-II, as can be seen from the relatively 
high value of Ax^ in table |T1 



3. Separating accretion and cooling-improvements of 
Sect. UTTbI 

In Sect. IIIIBI we discussed a procedure ('time shift') 
to separate temporally the accretion and cooling phases. 
The time shift produces two families of best fit values: 
when the 1MB data fall in the accretion phase (ij^^ ~ 
0) the positron temperature Ta has to increase; when 
the 1MB data fall in the cooling phase [tf^-^ ~ Tq) the 
temperature Ta can remain relatively low. We show in 
Fig. [5] the likelihood profile of Ta parameter, where the 
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above situation becomes evident.^ The figure refers to 
final model that includes also neutrino oscillations, but 
the same structure arises already as soon as the time-shift 
is included. The other local maximum of the likelihood 
besides the one shown in Tab. |T]is at: 

i?e = 10km, = 5.28MeV, = 4.74 s, ,^7) 
Ma = 1.27 Mq, Ta = 1.65 MeV, = 1.4 s 

with = 0.93 s and Ax^ = 10.2. The values ob- 
tained with this best-fit solution is very near to best-fit 
value shown in Tab.[Tl the difference being only ~ —0.5. 
Therefore this solution cannot be discarded on statistical 
basis. Thus, let us examine the physical content of this 
solution. When the data of 1MB belong to cooling phase, 
as in this solution, the value of Ta diminishes to account 
for the mean energy of the first KII events, and the initial 
accreting mass Ma increases to achieve the right number 
of detected events. This implies that the family of solu- 
tion with a temporal shift different from zero has larger 
values of Ma- The solution of Eq. [27]has a time constant 
Ta two-three times larger that the expectations and, more 
importantly, it has a value Ma twice the outer core mass; 
instead, the solution of Tab. U has a completely accept- 
able value of Ala- Since we expect that only a fraction of 
the outer core mass is exposed to the thermal positron 
flux, we are led to believe that the latter solution is more 
plausible than the former. 

4- Effects of oscillations-improvements of Sect. \III Q 

In the last line of Tab. [H we complete the parameteri- 
zation of the fiux by including neutrino oscillations. The 
solution is very similar to the one described in the previ- 
ous section and the astrophysical parameters are rather 
similar to the expectations. The total energy emitted by 
neutrinos in each emission phase is: Ec = 1.8 • 10^^ erg 
and Ea = 4.8 • 10^^ erg and the total binding energy is 
Eb ~ 2.2-10^'^ erg. The best-fit values have been obtained 
for normal hierarchy and with the assumptions discussed 
in Sect. IIII CI The flux of De that reaches the detectors is 
a combination of the and emitted within the star. 
The best-fit values showed in the table refers to the radius 
and temperature of De, that however are closely related 
to the i>f^ values. In fact the temperature of emission is 
TciPf,) = 5.5 MeV (due to Eq. ^ and the radius of z^^ 
neutrinosphere is Rc(v^) = 10 km (due to equipartition) . 

The assumption that we can neglect the during ac- 
cretion implies that the accretion flux is suppressed by 
the factor P — cos^ 612 ^ 0.7, so the best-fit of initial 




FIG. 2: Profiles likelihood for Ta parameter of accretion with 
the complete emission model (III.C). The red line is the family 
of solutions with t°f,g ~ 0.5 s. The dark line is the other family 
of solutions with t°^ = and Ma < IMq. The dotted line is 
the unphysical family of solutions discussed in Sect. B\ with 
Ta = 0.2 s. 



accreting mass has to increase by 1/P to maintain the 
same fit. This quantitative change has some impact on 
the interpretation of the multiple solutions: In fact, the 
family of solutions with t°^^j^ 7^ s and with lower values 
Ta (red line of Fig. H]) turns out to be characterized by 
values of Ma greater than 1 AIq. Thus, physical consid- 
erations on the meaning of the accreting mass suggest, 
as more plausible, the family of solutions shown with a 
dark line in Fig. [2] and in the last line of Tab. [J 

For inverted mass hierarchy, the effects oi — i' in- 
teractions are known to be relevant but their quantita- 
tive treatment is still under discussion, see e.g., [tsI. FtI]. 
This consideration, already, prevents us to draw firm 
conclusions. However, adopting the formulae given in 
Sect. IIII Cl to illustrate which are the possible effects, we 
can distinguish 2 main cases: large values of ^13, namely 
6*13 > 0.5° and small values of 613, namely 6*13 < 0.1°. 
In the first case, ~ and the survival probability 
P = U^i; thus, we find the same results as for normal hi- 
erarchy. In the second case, the flip probability is P/ ~ 1 
and P = J7g3 in Eq. The suppression of the Pe sur- 
vival probability P ~ 0, along with the assumption that 
the flux of and Dr are very small during accretion, 
implies the absence of electron antincutrino events dur- 
ing accretion^. Such a hypothesis can be tested with- 



* In technical terms, the case when there are multiple maxima of 
similar quality is termed as pathological likelihood. The statis- 
tical concept of 'pathological' solution should be distinguished 
from the concept of 'unphysical' solution; e.g., the one corre- 
sponding to dotted line of Fig. [2] and explained in Sect. IIII Bl 



Note that neglecting v — v interaction, the condition P ~ turns 
out to be realized in the other case, namely for large values of 

ei3 [zi. 
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and, to some extent, excluded by-SN1987A observations 
[79| . All in all, the result for inverted hierarchy should 
be taken with great caution for it depends crucially not 
only on the present understanding of the effects of v — v 
interactions but also on the incomplete description of the 
flux of Dfj, and Vt during accretion. 



5. Remarks on the two component models 

Having calculated the best fit values of the two com- 
ponent model, and having understood their meaning, we 
pass to discuss: (a) the evidence for the accretion phase 
(namely, we compare the two component and the one 
component model); (b) the errors on the best fit param- 
eters for the improved model; (c) the difference between 
the best-fit ECTA model and our improved model. 

a. Evidence for the phase of accretion In order to 
test whether the Ax^ of the models with accretion 
(Tab. m are just an effect of fluctuations, we perform 
a standard likelihood ratio significance test^ [SO]- For 
the four models of table HI wc find that we can reject the 
null hypothesis (=no accretion) with a significance level 
of a = 0.2%, a = 1.1%, and a = 2.0% (for the last two 
models) respectively Three remarks are in order: 

(1) The result a = 0.2% basically agrees with the claim 
of Lamb and Loredo: the ECTA model for emission, if 
correct, would lead to an important evidence for accre- 
tion. Our improvements in the likelihood (described in 
Sect, ini) do not change this inference significantly. 

(2) Also the other models permit to exclude the 'null 
hypothesis' that we test (namely, the absence of an ac- 
cretion phase) with the conventional 5% criterion. But 
the evidence becomes a bit weaker and if one prefers to 
be very conservative, this could suggest caution. This 
outcome can be easily understood by the fact that our 
model for accretion is more constrained and can account 
for certain features of the data (such as the difference of 
1MB and KII energies) only at the price of some tension, 
that is reflected by the increased value of a. 

(3) Obviously, even the conservative attitude does not 
forbid us to use the SN1987A data to learn something 
on accretion. It is the question that we formulate that 
changes: if we assume that the accretion phase exists, we 
can ask the data to determine the model parameters. 

b. Errors on the parameters The la errors obtained 
by a conventional, Ax^ — 1, Gaussian procedure [13] are: 

Rc - mtl km, Ma = 0.22t° j« Mq, 

Te = 4.6l°:^ MeV, = 2.4t°-j MeV, (28) 

r, - 4.71};^ s, r, - 0.55^0:^8 s. 




0.5- 
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Tc(MeV) 

FIG. 3; Two dimensional confidence regions for cooling pa- 
rameters Rc and Tc of with the complete emission model 

mm . 

The 1 sided, la errors for the offset times (obtained by 
integrating the normalized likelihood profile [13]) are:^ 

= 0.+°-°^ s, = 0.+"-^<' s, = 0+"-^^ s. (29) 

The couples of parameters that are more tightly corre- 
lated between them are Tc with i?c, and Ma with T^. 
In Fig. [3] and Fig. 2] we report the two dimensional confi- 
dence regions for these couples of parameters showing the 
90% and 68% contour levels; the correlations are quite 
evident from the figures. In these figures we focused on 
the family of solutions with = after testing that the 
other maxima have a very large value of accreting mass, 
as discussed in Sect. IIV B 41 

c. Stability of the best- fit values In order to show the 
stability of our result we investigated: 

1) different values for the exponent m = 1, 3, 4 (rather 
than m = 2) in Eq. [TH that describes the temporal be- 
havior of positron temperature during accretion; 

2) the value fc = 10 (rather than fc = 2) in Eq. [16] that 
describes the sharpness of the transition between accre- 
tion and cooling phases; 

3) deviations from the hypothesis of equipartition, by in- 
creasing or decreasing the ratio between the v^. and Dx 
luminosities by a factor of 2 |78j • 

In all cases, the changes less than 1 with respect to 
the best fit result; furthermore, the best fit values of the 



We calculate a = Y^i-i / exp(— a;^/2)/\/27r(ixi for X/T-i ^ 
Ax^ and taking Ax^ from Tab. U where the number of random 
variables Xi equals the new degrees of freedom v = 



^ Only 1MB had a reliable measurement of the absolute times; 
thus, we can use the time of its first event along with tf^Jg to 
infer the moment of the beginning of the collapse, presumably 
coincident with the emission of an intense gravitational wave. 
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FIG. 4: Two dimensional confidence regions for accretion pa- 
rameters Ma and Ta of with the complete emission model 
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Er [MeV] 




\mc\ 


Kl 


= 0.0 


20.0 


0.69 


1.00 


K2 


0.107 


13.5 


0.92 


1.00 


K3 


0.303 


7.5 


0.93 


0.81 


K4 


0.324 


9.2 


0.95 


0.93 


K5 


0.507 


12.8 


0.89 


0.84 


K6 


0.686 


6.3 


0.66 


0.10 


11 


= 0.0 


38 


0.02 


1.00 


12 


0.412 


37 


0.02 


0.55 


13 


0.650 


28 


0.05 


0.48 


Bl 


= 0.0 


12.0 


0.94 


0.99 


B2 


0.435 


17.9 


0.70 


0.85 



TABLE II: Accretion probabilities for the events occurred in 
the first second; all the other events have a probability to be 
due to accretion lower than 5%. The first three columns iden- 
tify the individual events. The last two columns are the prob- 
abilities that an event is due to accretion. The model used in 
the fourth column includes the improvements of Sect. and 
following LL sets Ma = 0.5 Mq; the one used in the fifth 
column includes also the improvements of Sect. CI 



astrophysical parameters change only within their 1 a 
errors given in Eq. 1281 

d. On the differences with the parameterization by 
Lamb and Loredo To iUustrate better the difference be- 
tween the ECTA model and our final emission model, we 
show in Tab, im the probabilities of the individual events 
to be due to accretion. A direct comparison with LL 
results [3^ is not possible, since a similar table is not 
given there; thus, we repeat their calculation following 
their prescriptions and make reference to the model LL* 
described in Sect. IIV B II The column LL* shows that 




<E>4MeV} 
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FIG. 5: i'e mean energy as a function of the time in the LL 
ECTA model (left panel) and in the \III model (right panel). 
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FIG. 6: i'e luminosity in the LL ECTA model (left panel) ana 
in the \IIld model (right panel). 



the early KII and Baksan events are due to accretion and 
those by 1MB to cooling. This proves that, assuming the 
Lamb and Loredo ECTA model, the fit takes advantage 
of the fact that the assumed energy distribution is 'com- 
posite': the low energy events of KII are explained by 
the accretion component, the high energy events of 1MB 
instead by the high energy tail due to the cooling com- 
ponent. 

The results of our model are shown in the column de- 
noted by nil CI Since our model has (by construction) a 
quasi-thermal spectrum at any time, both KII and 1MB 
early events have a large probability to be due to accre- 
tion. The tension between the diff'erent energies of KII 
and 1MB events leads to a slightly worse (see Tab.|I|. 
Finally, we compare in Fig. [5] the mean energy of De in 
the LL* accretion model, left panel, with the trend of the 
mean energy in our model [Til CI right panel, whereas in 
Fig. [HI we plot the luminosity obtained with the LL* 
model (left panel) and the model of Sect. IIII CI (right 
panel). The features of these curves are similar to those 
found in typical numerical simulations, see e.g. 113], with 
the exception of the average energy curve of the LL model 
that has a very pronounced jump. 



V. SUMMARY 

We presented an improved analysis of observations of 
SN1987A by Kamiokande-II, 1MB and Baksan. We recaU 
the main points: 

1. We collected in Sect. |lT]a large number of technical 
improvements: new detection cross section, procedure to 
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include the information on the direction of the events, 
treatments of the background and of the efficiency. The 
most relevant improvement is the last one, but none of 
these changes affect SN1987A data analysis in a crucial 
manner. 

2. We described in Sect. lIII Cl the various effects of neu- 
trino oscillations: the oscillations in the star (including 
the self interaction of neutrinos) and the Earth matter 
effect. We verified that oscillations with normal hierar- 
chy do not affect the results in an essential manner and 
discussed the inverted hierarchy case. 

3. We proposed in Sect. Ill II a new parameterization of 
the flux of electron antineutrinos emitted in the accretion 
phase. We improved on the energy spectrum, and most 
importantly, on the time-distribution: (a) We prescribed 
the temperature of the positrons to increase in such a 
manner that at i ~ Tq the average energy antineutrinos 
is approximatively continuous, that overcomes a short- 
coming of the parameterization of [s^ noted in and 
recalled in the Introduction, (b) We also prescribed that 
the number of neutrons exposed to the positron flux de- 
creases in time more smoothly than as in 36] ; in this way 
the luminosity is also continuous, as expected on general 
basis, (c) Finally, we avoided the simultaneous presence 
of cooling and accretion antineutrinos by time-shifting 
(delaying) the cooling phase of an amount t ^ Ta, again 
improving on [36| . 

4. We demonstrated that the improvements on the pa- 
rameterization are the most important. The most strik- 
ing feature is that the best fit parameters are rather sim- 
ilar to the expectations of the Bethc and Wilson sce- 
nario. This is in contrast with what happens using the 
Lamb and Loredo parameterizations of the flux, where it 
is needed to impose a prior in the analysis to avoid best-fit 
values outside the physical ranges. Furthermore, our flux 
leads to smooth luminosity curves and average energies. 
[The most appealing outcomes of their analysis, such as 
the neutrino sphere radius resembling the neutron star 
mass, the duration of the cooling phase ~ 0.5 s, the total 
amount of emitted energy, are practically unchanged.] 

5. We evaluated the errors and the correlations on the 
parameters and we can rule out the hypothesis of a one- 
phase model with a significance of 2%. 

From these calculations one draws two main mes- 
sages: an accurate choice of the model for the analy- 
sis of SN1987A observations is important; the agreement 
between the observations and the conventional expecta- 
tions is more than encouraging in the proposed model. 
We hope that these results will help to progress further 
in the understanding of this epochal observation. 
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APPENDIX: ACCRETION ENERGY SPECTRUM 

We derive the equations we used to describe the energy 
spectrum of accretion Vf.. 



1. Derivation of Eq. 1111 

We assume that, during accretion, the antineutrino 
flux is mostly produced by the weak reaction e'^n — > VeP- 
The rate of this reaction is given by: 



(3eC / dE, 



dE, 



(A.l) 



where iV„ is the number of target neutrons, assumed to 
be at rest, and where f3e is the positron velocity in natural 
units. The second integral yields the cross section as a 
function of the positron energy Ef,+ . The distribution of 
the positrons (ine+ — 2d^pp/h^ /(I +ex.p \(E^+ — Mp+ ) /Tg] ) 
has a negligible chemical potential [ij], thus 



dn^+ 

JeZ 



{hcf 



9e+{E^+,Ta) 



(A.2) 



where g^^ is given by Eq. 1121 The range of integration of 
E, in Eq. ED can be easily found from the expression: 



1 + {Ee -Pe C0S(/))/m„' 



(A.3) 



where cos (p = fie ■ is the cosine of the angle be- 
tween the directions of the positron and the antineutrino, 
COS0G [— 1, 1], and where 5+ = (m^ — mp-|-mg)/(2m„) — 
1.293 MeV « S- in Eq. [H The dependence on the an- 
tineutrino energy E,^ is as follows: 



dV 

dE, 



Sttc 

37^ = TTT^^n / dEe+ I3i 5e+ 



{hcf"" J dE, 
The cosine can take all values for fixed E, > E„ 
nie + (5+ 



E„ 



1 -t- nie/mr. 



1.803 MeV, 



(A.4) 
with 
(A.5) 



[for smaller E, , only certain values of cos (f) around 
cos = — 1 are allowed; but this happens in an inter- 
val of E, wide only about 1 eV]. The range for Ee+ is: 



Ee+ = i+A ' vA.6j 

where e — E,/mn and 1 -f A = (1 — e)^ — cos 0^. From 
this equation it is clear that the Ee+ range is pretty nar- 
row for the energies E, <^ rUn in which we are interested; 
thus, we approximate the expression in Eq. IA.4I as: 



« llfF^"5e+ {Ee+ {E,l To) JdEe^Pl 



(A.7) 
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where the positron distribution is calculated at the cen- 
tral point of the interval of cosine, namely at cos </> = 0: 



1 - E^/m„ ' 



(A., 



Eq. IA.7I gives the Pe flux in Eq. [11] the integral will be 
discussed a while. The advantage of Eq. lA.Tl is that the 
dependence on the parameter Ta has been extracted from 
the integral. The integral can be calculated once forever 
and we are left with a simpler expression. We checked 
that, in the most relevant range Ta = 1 — 4 MeV, the 
approximated expression agrees with the correct one at 
the 1% level for energies E^, < 10 Ta, and even better 
when dT/dEi, is integrated in E^: indeed, the rate F is 
precise at 0.1%. 



2. Derivation of Eq. [13] 

The cross section in Eq. I13[ precisely defined as a nu- 
merical approximation of the integral in Eq. IA.7[ was 
obtained adapting the calculation of [56] for the inverse 
beta decay reaction. In fact, the differential cross section 



[Ee 



ZOOTT nin Pi 



(A.9) 



has the same matrix element |A1p(s — u,t). What 
changes is the expression of the invariants, now given by: 



2mn{E^ + Ee) + ml and t 



-m;^ + 2nin{E^' 



Ee). The factor r{Ee) describes the small QED radiative 
corrections; we use expression in [56j . The constant in 
front to the differential cross section is 2 times smaller 
than the one for inverse beta decay, because the antineu- 
trino has 1 hclicity state whereas the positron has 2. The 



characteristic l//3e behavior of an exothermic reaction 
(such as e'^n — > D^p) is compensated by the 2 explicit 
factors (3e from the positron phase space and from the 
relative velocity between e"*" and n in the reaction rate, 
included in (Te+n- 

Our cross section compares well with the approxima- 
tion of Tubbs and Schramm l8l|, HI: 



TS 



= 1.7 X 10" 



35i 



E, 



(A.IO) 



with gA — —1.27, since the percentage deviation 100(1 — 
crj^/c^e+n) at E^ = 5, 10, 20, 30 MeV is just -1%, -2%, 
-6%, -11% (or -1%, -3%, -7%, -10% when comparing with 
the approximation of Ue+n in Eq. [THjl . 

The cross section used in [s^ is formally less correct, 
since it is the same as the above approximation but re- 
placing g\ —I- \gA\ ~ 1.254. [This is stated in Eq. (4.5) 
of LL and can be checked by the value of the energy ra- 
diated during accretion, their Eq. (6.2)]. The deviation 
100(1 - <y^+^l(Je+n) is not large; for E^ = 5,10,20,30 
MeV is 17%, 16%, 13%, 10% (or 18%, 16%, 13%, 10% 
when comparing with Eq. I13p . 

Our parametrization of D^. spectrum, Eq. I12( differs 
also for another reason with the one of LL, since we 
adopt the positron flux calculated in Ef.{E^) (defined 
in Eq. IA.8|1 . whereas LL use the antineutrino flux cal- 
culated in E^, namely, g El/[1 + exp(i^j^/Ta)]. Also 
this modification acts in the direction of increasing the 
expected flux. The difference can be quantified by eval- 
uating the integral of the fluxes $p^(Tq) = J ^^^^ dE^,: 
indeed, 1 - $^/'(T„)/$p^(Ta) = 54%, 35%, 26% or 19% 
for = 1, 2, 3 or 4 MeV; namely, our flux is significantly 
larger. 
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